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Abstract 

We study the Higgs boson effects on third-generation squark-pair production in proton-proton 
collision at the CERN Large Hadron Collider (LHC), including tt*, tb*, and bb*. We found that 
substantial enhancement can be obtained through s-channel exchanges of Higgs bosons at large 
tan/3, at which the enhancement mainly comes from bb, be, and cb initial states. We compute 
the complete set of electroweak (EW) contributions to all production channels. This completes 
previous computations in the literature. We found that the EW contributions can be significant 
and can reach up to 25% in more general scenarios and at the resonance of the heavy Higgs boson. 
The size of Higgs enhancement is comparable or even higher than the PDF uncertainties and so 
must be included in any reliable analysis. A full analytical computation of all the EW contributions 
is presented. 

PACS numbers: 12.60.Jv, 14.80.Da, 14.80.Ly 
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I. INTRODUCTION 



Supersymmetry (SUSY) is one of the most promising extensions of the standard model 
(SM). Not only does it provide a natural solution to the gauge hierarchy problem, but also 
gives a dynamical mechanism for electroweak symmetry breaking and a natural candidate for 
the dark matter. The simplest and most popular realization of supersymmetry is the minimal 
supersymmetric standard model (MSSM) [H, 0, 0]- The MSSM predicts the existence of 
scalar partners to all known quarks and leptons, via which the electroweak scale is stabilized. 
Since none of these SUSY partners have been found, the SUSY must be broken in our present 
world. Naively, we expect the SUSY particles to be heavier than their SM counter parts; 
however, naturalness arguments suggest that the scale of SUSY breaking, and hence the 
masses of the SUSY particles should not exceed C(lTeV). 

There exist some lower mass limits on these scalar SUSY partners. One of them comes 
from the direct search at the Tevatron. Current lower limits on the first- and second- squark 
masses are 200 — 300 GeV, depending on the gluino mass and the neutralino mass. An 
indirect limit comes from the Higgs mass bound of 114.4 GeV jij. It is well-known in the 
MSSM that radiative corrections can lift the Higgs mass at its tree-level bound (mz) to the 
current mass bound or more 0. The major correction comes from the top-stop loop. The 
current Higgs mass bound demands either (i) the top squark mass to be of order 1 TeV, or 
(ii) the mixing between the left- and right-handed top squark to be strong. While the first 
condition makes the search at the LHC experiments very difficult, the second option becomes 
very interesting. Not only can it satisfy the Higgs mass bound, but also allow a relatively 
light top squark as light as 200 GeV, which can certainly be produced at the LHC or may 
be even at the Tevatron. The left-right mixing effect in the first two generations of squark 
is negligible. Furthermore, the third generation of scalar fermions, t, b, and r are expected 
to be lighter than the corresponding scalar fermions of the first and second generations in 
Grand Unified SUSY models, because of the large Yukawa-coupling evolution. Therefore, 
potentially the top and bottom squarks are among the first SUSY particles to be discovered 
at the LHC. In this work, we focus on the third generation squark-pair production. 

There have been many works on hadronic production of top and bottom squarks [6[ . The 
leading order (LO) cross section for diagonal top-squark pairproduction (i=l,2) via gg 
scattering in hadron collisions was first calculated in Ref. [7|. However, these scattering 
reactions only lead to squark and antisquark pairs of the same flavor and same mass eigen- 
states, i.e., tit*, i 2 t*,, not only because gluons do not couple to t L i* R pair, but also because the 
coupling strengths to LL and RR are the same. Furthermore, such squark-pair production 
via QCD is weakly dependent of tan /3. The next-to-leading order (NLO) calculations have 
been improved by a number of authors including SUSY-QCD corrections. Recently, 

NLO electroweak contribution of 0(a 2 s a) as well as the lowest-order 0(a s a + a 2 ) electroweak 
terms to the production of diagonal squarks were studied and shown to be sizable [Io|. It 
was shown in Ref that the tree level electroweak contributions to the production of 
squark pairs at hadron colliders, which includes s-channel gauge boson exchanges as well as 
t and/or u channel gaugino exchanges, are comparable to the dominant QCD contributions 
in some cases. Note that in that study [ll[ the bb initial-state contribution was not included. 

Non-diagonal squark-pair production, like tit*, and bib*,, is possible at tree level via Z 
boson exchange, as the LL and RR couplings are different Mixed top and bottom 



squark production tib* is also possible via an intermediate W boson in the 2 — > 2 subprocess 
or via gW fusion in the 2 — > 3 subprocess [l2j|. All these non-diagonal trf*, tib*, and bib* 
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pairs with % ^ j are produced via exchanges of electroweak gauge bosons, and therefore the 
production rates are smaller than the corresponding diagonal pairs. 

In this work, we point out that the s-channel exchange of various Higgs bosons plays an 
important role here, namely, it can substantially enhance the production rate, especially at 
large tan (3. In this paper, we investigate the importance of electroweak channels including 
the neutral Higgs bosons for trf* and b^b* pair production and charged Higgs boson for 

tib* pair at hadron colliders when one or both of the initial-state partons are the bottom 
quark. We anticipate the contributions from the Higgs exchanges will be substantial at large 
tan (3 region, where the smallness of the bottom-parton luminosity can be compensated by 
enhancement of the bottom Yukawa coupling. 

Furthermore, it also allows resonant Higgs production for relatively heavy Higgs bosons 
H ± , H°, and A . In fact, with m^o^o #± > + m^, the nondiagonal squark pair pro- 
duction offers an interesting possibility to study the squark-squark-Higgs couplings right at 
the Higgs boson resonances. In addition, with the intermediate charged Higgs boson all left- 
and right-handed squark pairs t^^b^R can be produced while only t L b L can be produced via 
the intermediate W exchange. Thus, we anticipate the production via intermediate charged 
Higgs boson could be dominant in some region of parameter space. We explore the MSSM 
parameter space relevant for our study. 

The organization is as follows. In the next section, we will write the details of the cou- 
plings and mass matrices. In Sec. Ill, we describe the formulas for squark-pair production. 
In Sec. IV, we discuss the effect of phenomenological constraints on SUSY parameters, 
followed by numerical results on production cross sections. We conclude in Sec. V. 



II. FORMALISM 



We start with the following superpotential 

w = e ab >;y/;//:;r; - v %Q)u h d Di - vI^h^e* + ^n*] , (i) 

where ei2 = — 621 = 1, i, j are family indices, and y u and y d represent the Yukawa matrices for 
the up-type and down-type quarks, respectively. Here Q, L, U c , D c , E c , H u , and Ha denote 
the quark doublet, lepton doublet, up- type quark singlet, down-type quark singlet, lepton 
singlet, up- type Higgs doublet, and down- type Higgs doublet superfields, respectively. 



A. Quark mass matrices 



The Higgs doublets develop vacuum expectation values (VEV), which break the elec- 
troweak symmetry, are 

where v ~ 246 GeV, and the quark mass terms are given by 

vsin/3 . v cos f3— 

£ quark = -Uij ^2 u m u Lj - Uij ~f= d Ri d Lj + h.c . (3) 
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The Yukawa matrices y u and y d are diagonalized by bi-unitary transformations: 














= Dl,r 1 




1 • 1 




= u LiR 


L,R 




L,R 




L,R 




(4) 



L,R 



and 



Ui L^)u L =( 7 11), D t Dl = (o m, ] . (5) 

V / \^ irii/ \ J \ m b 

We define the diagonal Yukawa matrices Y u and Y d as 

/^/m„00\ /9/m d 00 
1™ = -4^= m c , Y d = m s | . (6) 

usm ^ \ m t / wcos ^ \ m fe 

In order to avoid excessive flavor-changing neutral currents, we assume that the squark 
mass matrices are in alignment with the quark mass matrices, i.e., they are diagonalized by 
the same bi-unitary transformations: 










f d A 


= Ul,r 1 




1 • 1 




L,R 




L,R, 






= Dl,r 5 , (7) 

L,R \ b ' L ,R 

Here di,d 2 ,ds are in interaction basis while d,s,b are in mass eigenbasis. Without loss of 
generality, we can make the choice that the right-handed quarks and squarks are already in 
the mass eigenbasis, ie., Dr = Ur = I. The information on the left-handed unitary matrices 
Dl, Ul is encoded in the Kobayashi-Maskawa (KM) matrix as 

V KM = U ] L D L . (8) 



B. Squark mass matrices 

Even after we have rotated the quark and squark mass matrices into family- diagonal 
form with the same bi-unitary transformation, within one family the left-handed and right- 
handed squarks will mix due to the soft terms and F terms in the Lagrangian. This LR 
mixing is proportional to the quark mass concerned such that the LR mixing for the first 
two generations are negligible while it can be substantial for the third generation. From now 
on we only concern the LR squark mixing in the third generation. 

The squark mass-squared matrix in the L-R basis have the form 



with 



M- 



LL 

2* 



III 

ml 



LR 



LR " b RR 



(9) 



m 



LL 



711 



RR 



Ml + m 2 q + m| cos 2(3 ( If 
M\ + m 2 + m| cos 2f3 e q s 



2 

Wi 



(10) 
(11) 



ni 



LR 



m n 



A* - /i (tan/5)~ 2/ 3 



(12) 
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where J| = ±1/2 and e q are the third component of the weak isospin and the electric charge 
of the quark q. In Eqs. fl§]), /x is the supersymmetric Higgs mass parameter, M| the soft- 
breaking mass parameter for the squark iso-doublet (g^, q'i), and M R are the soft-breaking 
mass parameters for the iso-singlets q R . They can be different for each generation, but 
for simplicity we will assume equal values for all generations Ml = M R = Msusy i n our 
numerical analysis. A q are the parameters of the soft-breaking scalar three-point interactions 
of top- and bottom-squarks with the Higgs fields. 

The hermitian matrix in Eq. (Q is diagonalized by a unitarity matrix R q , which rotates 
the current eigenstates, qi and q R , into the mass eigenstates q\ and §2 as follows, 

<?i \ _ uq ( Ql 




yielding the physical mass eigenvalues, with the convention m qi < rriq 2 , 



m l,2 = \ Wll + mlj, =F ^{ml L -m 2 RR y + A\m 2 LR \^ . (14) 



The mixing angle 9 q obeys the relation 



tan 20,= 2 2m '* 2 . (15) 



m LL ~ m RR 



Hence, for the case of the supersymmetric partners of the light fermions, L-R mixing can 
be neglected. However, mixing between top squarks can be sizable and allows one of the 
two mass eigenstates to be lighter than the top quark. Bottom-squark mixing can also be 
significant if tan (3 is large. 



C. Higgs and gauge bosons interactions with quarks and squarks 

Let Hk = (h°, H°, A , G°) (k=l..A), one can write the relevant Lagrangian density in the 
(gi,g 2 ) basis as following form (2, j=l,2) 

^relevant = 9H k qqHkqq + (G k )ijH k q*qi + (G 5 )ijH + q*qi (16) 

+ gq{A\ m P R + BfjxUi + <7?'(4A + KP L )xl% 
+ H + u (Y u cos (3P l + Y d sin (3P R ) V KM d 

+ G + u (Y u sin(3P L - Y d cos(3P R ) V KM d 
+ h.c 

The Feynman rules are ign k qq for k =1,2 and 75(7^^ for k= 3,4 where we list only Higgs 
bosons couplings to quark b, and 

gm b sina gm b cosa gm b 

ghbb = 9mb = —z gAbb = -7; tan p. (17) 

2m w cos p 2m w cos p Imy/ 
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The couplings Gk and G5 are given in the mass basis by 

G k = K*G k (K*) T , (k = l..A), G 5 = K i G 5 (K i ) T (18) 

where Gk and G5 are the couplings in the {ql,Qr) basis, and their explicit forms are 
shown in Appendix El Finally, A q m , B\ m (m=1...4), and £?,, K q u (1=1,2) used in Eq. ffTBl) are 
defined also in the end of Appendix [A] 

The W and Z bosons interactions with the quarks are given by 

C = -JLu^V KM P L dW; - gZqY(C qL P L + C qR P R )q + H.c. (19) 
V2 

with 

C qL , R = — L_ (II - eg sin 2 6 W ) (20) 

COS Uyy 

The W and Z bosons interactions with the squarks are given by 

C = -iggwMjWffi&di - iggz M] Z~q*d^~ qj + h.c. (21) 



with 



1 / cos 6u cos #j — cos 6u sin 6^ 



i.i 



a/2 V — sin 6u cos #j sin 6q sin 6 d 

and 

/or cos 2 On — e„ sin 2 6 W —hlii sin 20, 



cos 6 W V - 1 /f L sin 20,- J| L sin 2 8 - - e q sin 2 W J . . 1 ] 



D. Radiative corrections to the Yukawa couplings 



Note that both initial and final states have the bottom Yukawa Y b dependence. It is 
now well established that the coupling of the scalar bottom b to the up-type Higgs doublet 
induces a modification of the tree-level relation between the bottom quark mass and its 
Yukawa coupling [3, 20, 21, 22|. Those corrections are amplified at lar ge t an / 3. The 
modifications can be absorbed by redefining the bottom Yukawa coupling as m, 13, EH, Q 



Y b 



V2 



tan (3 



v cos (3 v cos (3 1 + Aft v 1 + A5 
where the second term is valid for large tan (3 and the SUSY-QCD corrections lead to 



(24) 



A,, 



2a. 
3tt 



t\2 



(Y 

H nig tan (3 /(m^ , , m~) + V A tan (3 1(m Tl , , /i) 



where m q is the gluino mass, and the function / is given by 



I(a, b, c) 



1 



(a 2 - b 2 )(b 2 - c 2 )(a 2 - c 2 ) 



1 2 
fc 2 



, 2 lnS + foVln^ + cVln^ 

a 2 ' 



(25) 



(26) 



(3 



In Ab we only keep the dominant contributions from the gluino-sbottom and charged- 
higgsino-stop loops because they are proportional to the strong coupling and to the top 
Yukawa coupling, respectively, while neglecting those that are proportional to the weak 
gauge coupling. Note that A b is evaluated at the scale of SUSY particles M susy , where the 
heavy particles in the loop decouple, whereas the bottom Yukawa coupling Y b (Q) at any 
scale Q is determined by the running 6-quark mass m b (Q) at the scale Q: 

The contributions to the bottom Yukawa couplings which are enhanced at large tan/3 



can be included to all orders by making the following replacements |19l . I2d . 1231 . 124 



1 — A b / (tan j3 tan a) 

9hbb —> Qhbb : : a (2'J 

J- + L±b 

1 + A b tana/tan/3 . . 

QHbb — > QHbb : ; A (28) 

1- A b /tan 2 /3 

9 Abb -> 9 Abb YTA ' ' 

Y b ^<g> i (30) 

v cos (3 1 + A b v ' 

We now have all the tools to compute the production cross section of sbottom and stop. 



III. SQUARK PAIR PRODUCTION 

In this section we discuss squark-pair production. Let us define our notation for the 
convenience of the following formulas. The momenta of the incoming quark q and anti- 
quark q, outgoing squark (fj and outgoing anti-squark q* are denoted by pi, p2, k\ and 
&2, respectively. We neglect the quark masses of the incoming partons. The Mandelstam 
variables are defined as follows 

8 = (P1+P2) 2 = (h + h) 2 

t = (pi - hf = (p 2 - k 2 f = 9 ' 2 gj - S - (1 - cos 9*) 

U = (ft - k 2 f = ( P2 - h) 2 = 9 ' 2 9J - I (1 +0COS0*) (31) 

where /5 = A 1 / 2 (l, m~,/s, m~ /s) and #* is the scattering angle in the center-of-mass frame 
of the partons. 

A. Hadronic production of gg, — > tit^bib^ 

The production of top and bottom squark pair proceeds via the following qq- and gg- 
initiated subprocesses which are depicted in Fig. [T](a) - (h). 

QQ,99 -^Utphb*: (32) 

where = (1, 2). 
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(i) (j) (k) 

FIG. 1: Tree-level Feynman diagrams for squarks pair production via qq, gg, and qq' annihilation. 
With (i,j) = (1,2), q = (it, d, c, s, b), S = h°, H°, A , G° and V = 7, Z, g. In diagram (b), the photon 
and gluon do not contribute to non-diagonal production. 



Note that the photon and gluon do not contribute to production of of qiq^, because 
of electromagnetic or color conservation. Consequently, at tree level the above reactions 
proceed only through s-channel Z boson and Higgs bosons. If the initial state is uu, dd, ss 
or cc, the Yukawa couplings are so small that we consider only the contribution from the Z 
exchange diagram. While in the case of bb, the cross section will directly determined by the 
size of Higgs coupling to a squarks Sqifo where S = h°, H° or A and Zq~iq 2 couplings, which 
are proportional to sin 29q of the squark-mixing angle 9g. therefore, these processes can be 
used to probe the mixing angle 9q. The analytic expressions for bb — > tit* and bb — > bib* are 
given in appendix iBl 

The gluon-gluon fusion into tit*, and bib*, only goes through loop diagrams, as shown in 
Fig. [21 These one loop contributions are of the order of either a 2 s or a s a. 

The hadronic inclusive cross section for qid* production in proton-proton collisions at a 

total hadronic center of mass energy \/~S can be written as [29[ 

" x dm r 1 d£PP 



/j- dC - f d£ pp 

dT -^VLo(qq -> q$j){rS) + J dr—^-d LO {gg -> q$){rS) (33) 

where r = (m|. + m-^/S, and the parton luminosity is 

f ^ dx T T 

I — [f a (x,ii F )f b (-,iJ lF ) + f b (x,fi F )f a (-,ii F )] (34) 

J T x 1 + d ab L x x 



dCZ f 1 dx 1 
dr 



where f a (x, fi F ) is parton distribution functions (PDF) for each type a in the proton carrying 
a fraction x of the proton momentum at scale ji F — + ms.. 



S 




S i S i 



o 9 Q 

s\ v\ 

(t) 



5 £i 



F 



g F 12 
(u) 



FIG. 2: One-loop Feynman diagrams for squark-pair production in gluon-gluon fusion. With i ^ j 
= (1,2), F = g, x°, X ± , S = h°, H°, A , G° and V = y,Z, g. 

B. Hadronic production of du — > ijK 

In this case, the Mandelstam variables are defined as in Eq. (l3Tj) with = and 
= At hadron colliders, the production mechanism proceeds via the conventional 
Drell-Yan prcoesses with the charged Higgs boson, charged Goldstone and charged gauge 
bosons, as depicted in Fig. [1] (i)-(k). The s-channel diagram with charged Higgs exchange 
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dominates when \/I is close to M^±. The analytic expressions for the cross sections are 
given in appendix iBl 



IV. NUMERICAL RESULTS 

A. Phenomenological constraints 

In this section, and before presenting our numerical results, we would like to list the 



phenomenological constraints included when determining the allowed parameter space [13, 
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The most stringent constraint generally arises from Ap SUSY which receives contribu- 
tions from both stop and sbottom. The extra contributions to the Ap SUSY parameter 



from the stop and sbottom sector [15l . Il6l ] should not exceed the current limit from 



precision measurements |17|: Ap SUSY < 10 c 



The soft SUSY-breaking parameters A q at the weak scale should not be too large in 
order to keep the radiative corrections to the Higgs masses under control. In particular 
the trilinear couplings of the third generation squarks A t ^, they will play a particu- 
larly important role in the MSSM squarks/ Higgs sectors. These parameters can be 
constrained in at least one way, besides the trivial requirement that it should not make 
the off-diagonal term of the squark mass matrices too large to generate too low masses 
for the squarks. should not be too large to avoid the occurrence of charge and 
color breaking (CCB) minima in the Higgs potential. To avoid such minima at tree 
level, At b have to satisfy the following tree level conditions 18| • 



A 2 < 3 (m| + m\ - 2m 2 + ~m| cos 2(3 + M 2 H2 + p 2 ) , (35) 
A\ < 3 (ml + ml - 2m 2 h + X -m\ cos 2(3 + M 2 Hi + p 2 ) . (36) 

With Mff = (m 2 A o+m 2 z ) sin 2 (3— l/2m| and Mj^+Mj^, = m 2 AQ . The above constraints 
depend on p and A tjb explicitely. 

Another constraint which has been imposed is the perturbativity of the bottom Yukawa 
coupling. Since radiative corrections to the bottom Yukawa coupling have been im- 
plemented as in Eq. (1301) . the Yukawa may blow up when SUSY parameters varies. 
Thus, we restrict Y b < (Air) 2 . 

We have imposed also all the experimental bounds on squark, chargino, and neutralino 
masses as well as Higgs boson masses [l7| . 



B. Numerical results 

In this section, we present the numerical result for inclusive production cross section 
of diagonal, non-diagonal and mixed squarks at the LHC with a proton-proton center-of- 
mass energy cross sections of 14 TeV. In our numerical calculations the following SM input 
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parameters were chosen 
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mt = 171.9 GeV , m w = 80.398 GeV , m z = 91.1878 GeV 

G F = 1.16637 x 1(T 5 GeV" 2 , V cb = 0.04, m b (m b ) = 4.25 GeV (37) 



The running QCD coupling a s was evaluated at the two-loop level |25| and the CTEQ6L 
PDFs MM 

were used to calculate the various cross sections. Moreover, in order to improve 



the perturbative calculations, 1-loop running masses m b (Q) were taken as following: 

m b (Q) = mp(Q) = mp(Q) (l + ^ (38) 

where includes the SM QCD corrections. 

First of all, we investigate the effect of varying the MSSM parameters for which the 
correction to A b term are expected to have a large impact. In the limit where the squark 
and gluino masses have approximately the same value, denoted by the common SUSY mass 
M SUSY , the Eq. (1251) simplifies to 

I( m susY,m SU SY,msusY) = -z-^T • ( 39 ) 

Zm SUSY 

Furthermore, if \i is of similar size, the first term in Eq. (125j) is dominant and reduces to 

A b fa sign(/i)— 1 tan/3 (40) 
on 

So, for large tan/3 this effect can be 0(1) and does not vanish for a heavy SUSY spectrum. 
The sign of \i is the decisive factor in determining whether the corrections will enhance or 
suppress the cross section for the processes of pp —> bb —> bib*, tit*, Ub*. We assume 
the universality of soft SUSY breaking trilinear couplings: A t = A b = A T = Aq. We 
parameterize the squark sector using the following input parameters: tan /3, fi, A Q and the 
gluino mass m§. The MSSM Higgs sector is parameterized by the mass of CP-odd m^o and 



tan/3 as well as by M S usy, A b<t and fi for higher order correction s |26l . |27j. All the MSSM 
Higgs masses and parameters are computed with FeynHiggs code |26j. 

We present the results of diagonal squark production of bib\ in Fig. [3] and tit{ in Fig. IU 
respectively. At the LHC, the diagonal pair production of tit* and b\b\ is dominated by 
gg fusion. However, it is noted that the to-initiated subprocess already surpass the light 
qq- initiated channels when tan/3 ^17, because of the enhancement from the Higgs boson 
couplings to the bottom quark as well as from the large coupling of Higgs to a pair of squarks. 

The cross sections a(bb, gg, QQ — > Mi and bib\) are plotted as a function of pseudoscalar 
Higgs mass m^o for a large tan/3 in the right panel and as a function of tan/3 at around 
the resonance m^o w 2rriq, where q = b,t, in the left panel. The size of these cross sections 
depend strongly on the Higgs mass and tan/3. We note that indeed the leading order cross 
section a(bb — > t\t* and bib*) increases like tan 2 /3 for large values of tan/3. To understand 
this we separate various contributions to the production rate. The bottom-induced Drell-Yan 
contribution proceeds through an s-channel photon or Z boson. The coupling to the photon 
is independent of tan/3 and the dependence of the Z couplings is small for large enough 
values of the pseudoscalar Higgs mass. A second set of s-channel production processes is 
made possible by the incoming bottom quarks and their finite Yukawa couplings to CP- 
even Higgs bosons h°, H°, but the CP-odd pseudoscalar exchange is forbidden by the CP 
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symmetry of the final state. The couplings h°bib\ and H%ib\ depend strongly on tan/? and 
on other MSSM parameters such as \x and soft trilinear terms A t ^. Note that the dependence 
on tan (3 and fi parameter also comes in through A& corrections. 

The relative size of the Higgs contribution compared to Drell-Yan process show up at 
large tan (3 and near the resonance. Note that the interference between Higgs contributions 
and the QCD are zero due to color structure. The process bb — > q\q\ receives EW contri- 
butions from t-channel exchange of a neutralino and gluino if q =b and chargino exchange 
contribution if q = t. The corresponding Feynman diagrams are shown in Fig. [Tfc) and (d). 
Combining s-channel and t-channel including Higgs contributions, we see that they interfere 
constructively. Since QCD contributions dominate even after inclusion of the electroweak 
diagrams, the overall behavior of the total cross sections does not change ulucIl For exam- 
ple, for tan/3 ~ 20, the EW increase the cross section for the pair production of b\b\ squarks 
by about 2 pb, whereas they only contribute 24% to the total cross section. 

In gluon-gluon fusion, the non-diagonal squark ^g* with i ^ j cannot be produced at the 
lowest order but via loop diagrams (see Fig. 2). The production cross section is therefore 
of the order 0(a^ + a 2 s a 2 ). This higher-order cross section will be small compared to qq 
tree level contribution. The 66-fusion into non-diagonal squark pair proceeds via CP-even 
Higgs bosons h°, H°, and CP-odd Higgs boson A in s-channel, gluino and neutralino in 
t-channel for b\b 2 and chargino in t-channel for t\t 2 pair production. The cross sections 
<t(&i&2 + b\b X) and a(tii* 2 + t 2 t{) are plotted as a function of tan/3 in the left panel and 
as a function of in the right panel of Fig. [5] and Fig. 0, respectively. As we can see, 
the cross sections for both processes are small due to the absence of tree-level gg fusion 
and suppression from the phase space. For intermediate values of tan/3 due to negative 
interference term between s-channel and t-channel diagrams all processes are comparable in 
size. For large tan /3, A& effect induces a large enhancement for bb — > qiq 2 + <72<7i with q = b,t 
by an order of magnitude. The size of the to-initiated subprocess is already above the gg or 
gg-initiated channels for tan/3 > 5, as shown in Fig. [5] and Fig. [61 Note that in the case of 
bb^ b\b 2 , we can see a spectacular enhacement for large tan/3 (see Fig. left). The origin 
of this enhancement is due to a factor of tan/3 in the bb A coupling and another factor of 
tan (3 in the A°bib 2 coupling, which results in tan 4 (3 enhancement in cross section. This is 
a very interesting result that we cannot ignore the to-initiated subprocess with s-channel 
Higgs boson exchange. 

Table [D shows the cross sections including QCD, EW and Higgs effects. It is clear that 
Higgs effects enhance the cross sections for negative \x by one to two orders of magnitudes. 
While they reduce the cross sections for a positive \i. Hence, we can see that the EW and 
QCD contributions are comparable for diagonal pair production 

Stop-sbottom pair production is dominated by the W exchange diagram with ud-initiated 
subprocess. The charged-Higgs contribution through to-initiated subprocess can get to com- 
parable size of cross section at very large tan (3 > 35, as shown in Fig. Note that the 
charged-Higgs couples to the cb quarks with an enhancement from tan/3 but also with a 
suppression from Vcb- As we can see, large tan (3 limit overcome easily the V c b suppression. 

V. CONCLUSIONS 

We have reviewed the theoretical status of squark pair production at the LHC. The 
evaluation of the full electroweak contributions has been described in details. The Higgs 
effect enhance the cross section by about 10-25%, which is comparable to the size of NLO 
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tan/3 niA® (GeV) 

FIG. 3: Production rates for b\b\ pair production as a function of (left) tan/3 and (right) m^o. 
Other SUSY parameters are chosen to be M SU sy = 490 GeV, M 2 = 200 GeV, m~ a = -[i = \ TeV, 
A) = H40 GeV. 
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FIG. 4: Production rates for t\t\ pair production as a function of (left) tan/3 and (right) m^o. 
Other SUSY parameters are chosen to be M SU sy = 350 GeV, M 2 = 200 GeV, m 5 = 550 GeV, 
H = -960 GeV, A = 750 GeV. 
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FIG. 5: The cross sections of pp — > bib?} + b%b\ pair production as a function of (left) tan/3 
and (right) m A o. The SUSY parameters are chosen to be Msusy = 400 GeV, Mi = 200 GeV, 
m-g = 1000 GeV, p = -500 GeV, A = 500 GeV. 
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FIG. 6: The cross sections of pp — > tit^ + ^1 pair production as a function of (left) tan/3 and 
(right) m^o . The SUSY parameters are chosen to be M S usy = 400 GeV, M 2 = 200 GeV, m~ g = 1 
TeV, // = -500 GeV, A = 500 GeV. 
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0.87 
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TABLE I: The effect of the sign of the \x parameter on the production cross sections (in fb). The 



SUSY parameters are chosen to be Aq 



1 TeV, Msusy = 490 GeV and M 2 = 200 



GeV, tan (3 = 20 and the Higgs masses are taken at the resonance. 
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FIG. 7: The cross sections of mixed squark pp — > pairs production as a function of (left) tan/? 
and (right) m H ±. The SUSY parameters are chosen to be Msusy = 200 GeV, M<± = 200 GeV, 



1000 GeV, = -200 GeV, A = 200 GeV. Note that for non-diagonal production qq> 



we have taken into account their hermetic conjugate. 



fb* 



correction as well as to the PDF uncertainties. In the case of non-diagonal squarks pro- 
duction, we have seen some enhancement for large tan/3. Those processes can be used to 
extract some information on the squarks mixing angles. One concludes that LO electroweak 
contribution has to be taken into account for any reliable prediction. 
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APPENDIX A: FEYNMAN RULES 



In this appendix, we give the relevant Feynman rules for our study. 



1. Higgs-squark couplings 

In this appendix we collect the couplings of the h°, H° and ^4° bosons to the squarks q~i 
with q= t, b and i= 1,2, which are relevant for our analysis, by using the mixing matrix 
which rotates the left- and right-handed squark fields q^ and q^, into mass eigenstates qi as 
Squark couplings to Higgs boson are given by where a is the mixing angle in the neutral 
Higgs sector, we have used c a = cos a, s a = sin a, c a +p = cos(a + (3), s a+ p = sin(a + /?), 
S/3 = sin/3 and cp = cos/3 as abbreviations. We note the following properties of the above 
couplings ghPfth = ghPMxi 9H^q 2 = QWMi and 9aom2 = -9A°q a q v Squark couplings to 
quarks and either charginos x or neutralinos x° are straightforward, but somewhat more 
complicated by the mixing angles associated with the x and x° niass eigenstates. 



a. Squark - Squark - h 



9 m z r c 
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(Al) 



for \ up > type squarks respectively, a is the mixing angle in the CP even neutral Higgs 

I down J 

boson sector. Y q are the Yukawa couplings: 

9 m t . g m b 



Y 



m w sin (3 



Y 



m w cos (3 



(A2) 
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b. S quark - S quark - H 
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Notice that G 2 can be obtained from G\ by the replacement a — > a + 7r/2, i.e 
and s a Cq . 

c. Squark - Squark - A 
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2. quark-Squark-neutralino and quark-Squark-chargino 

a. quark - squark - neutralino 



A g 1 / -4m w s p B% i N^ 1 + 3c w m n N^B$ L \ (m 

ZV2c w s p s w \ 2c p M w s w N* ml Rf 2 + 3c w m d N^Ri J ' 

g 9 = 1 f 3c w m u ^4-Rf 2 + m w spR^(s w N ml + 3c w N m2 ) 1 ^ Ag s 

Sy/^cwSpSw \ 3c w m d Rf 2 N m3 + c^m^i^f^s^A^i - 3c w N m2 ) J 

Here iV is the 4x4 unitary matrix diagonalizing the neutral gaugino-higgsino mass matrix [jj • 



b. quark - squark - chargino 



, J* I - f (2mws Rt 1 V kl - s/2m u R%V k2 ) . 
" " < ^ / ■ K " - - \ ^ { 2m w cM^ - V»n d RiU t2 ) > ' (A1 °» 

Here U and 1/ are the 2x2 unitary matrices diagonalizing the charged gaugino-higgsino 
mass matrix ll|. 



APPENDIX B: PRODUCTION RATES 

The production of squark pairs, as initiated by bb annihilation, involves gluon, photon, 
Z, W and Higgs bosons in the s-channel as well as gluino and neutralino exchanges in 
the t-channel. Since left-and right-squarks generally have different masses we present the 
differential cross section for each subprocess separately in the mass basis. The spin and color 
averages are taken into account. i,j = 1...2, m,n = 1...4 and k,l —1...2. 
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The differential cross section for b b — » bi b* (diagonal) is given by 
dfr(bb^bib$) 2tt 
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Where D^ 1 = s — m\ + im<$,T<$,. The imaginary part in the D^ 1 is the Breit-Wigner 
prescription for regulating the $ pole. 
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while the differential cross section for off-diagonal (27^ j) is given by 
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where summations over 
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a. bb — > U t\ 



Similarly, the differential cross sections for tit* and t(t* are respectively given by 

2a 2 



dt 9s 2 1 1 



+ 



e\a 2 a 2 2 



+ 



+ 



9sa' 



4(ttt — m~ ) 
9a 2 



9huu9hbbDh + g H t.j.gHbbD H 



i(t-m 2 )(t-m 2 ) 

Xk X-l 



2m~+m~+s 

(Afc) + + 7J7 ^ 

(ru — ra- 
ti ' 



-C£ik) 2 (£< 



— + -rtI^ + c bR)9zutA s - M z) D 



- 3a' 



(t-ml+) 

£ (4g|>g 

2 (t-m 2 ^) (tu-rrd) 

Xk *« 



6 s 

(t -u) fe t 3 



(B3) 



x 



9huu9hbb{s - M 2 )L> 2 + g HlA g Hbb (s - M 2 H )D 



3a 2 (s - Mf) ^-y, , r , r Fyn n ' ^ 2 [(4 fc ) 2 + (4) 2 ] 



2 

Xk 



2s (t - ml 



da(bb^tit*) _ 2tt 
+ 



a 



(B4) 



9sa^ 



4(t , u — m~m 2 
9a 



2 



Kt-m| + )(t-m| + ) 



QhtajQhbbDh + gmj )9HbbD H + g A t.tgAbbD A 



(4,) 4 + (4) 4 + tt: Xfc (4) 2 (4 V 



— 3a 



f + 3 



(i-u) 



g z uiC bL + C bR )(s-M z )D z 



2 (t — m 2 + ) 2 (tu — m~ m~ ) 
x - M 2 )D 2 h + grfftgmbis - M 2 H )D 2 H + ^.^ 66 (s - MD^J 



3a 2 (s - M 2 ) 
2 {t-ml + 

Xk ' 



[{C bL {K] k f + C bR (4 k ) 2 ]D 



21 



2. stop-sbottom production 
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